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Abstract 

We study a many-body mixture of an equal number of bosons and two-component fermions with a 
strong contact attraction. In this system bosons and fermions can be paired into composite fermions. 
We construct a large N extension where both bosons and fermions have the extra large TV degrees 
of freedom and the boson-fermion interaction is extended to a four-point contact interaction which 
is invariant under the 0{N) group transformation, so that the composite fermions become singlet 
in terms of the O(A^) group. It is shown that such O(A^) singlet fields have controllable quantum 
fluctuations suppressed by 1/7V factors and yield a systematic 1 /A/"- expansion in terms of composite 
fermions. We derive an effective action described by composite fermions up to the next-to-leading- 
order terms in the large TV expansion, and show that there can be the BCS superfluidity of composite 
fermions at sufficiently low temperatures. 
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1 Introduction 



The study of boson-fermion mixtures has a long history originating from the analysis of dilute solu- 
tions of ^He atoms in superfluid ^He 1 . For a weakly-coupled boson-fermion mixture, it is known 
that the density fluctuation of the bosonic background induces an attraction between the fermions, 
which enhances the transition temperature to the BCS superfluidity or leads to fermionic superfluidity 
even without a bare attractive potential between fermions |2,3 . On the other hand, in the strong 
coupling regime, it is possible to form bound states between bosons and fermions, called composite 
fermions (CFs) or simply dimers (tightly-bound molecules) 4 . Therefore, phase structures in the 
strong coupling regime may differ from those in the weak coupling regime, and it is expected that 
there occurs the superfluidity of CFs at low temperatures [5], which greatly motivates us to model 
superfluid hadronic matters in dense QCD in terms of boson-fermion mixtures where small size di- 
quarks correspond to the bosons, unpaired quarks to the fermions, and the extended nucleons are 
regarded as the CFs [6||8]. 

Recent developments in atomic experiments have made it possible to realize boson-fermion mixed 
gases in the laboratory. Atomic interaction between different species can be tuned with the use of 
Feshbach resonance techniques [£-11 . Recently, the formation of heteronuclear Feshbach molecules 
has been observed in a boson-fermion mixture of ^^Rb and ^^K atomic vapors in a 3D optical lattice 12] 



and in an optical dipole trap 13 



From a theoretical point of view, there are several non-perturbative studies on nonrelativistic 
atomic gases. In particular, the large N method provides a systematic expansion with the cor- 
responding diagrammatic representations, and the applicability of its results to the physical cases 
at = 1 can in principle be tested by systematic estimates of higher-order contributions. The 
transition temperature of the dilute interacting Bose gas has been calculated with the use of 1/N- 
expansion [14[|15|. Also, the l/A'-expansion for the nonrelativistic Fermi gases has been developed in 



Refs. 16 -18 . A review of large A^ expansions in O(A^) and U(A^) quantum field theories, which deals 
with non-perturbative aspects of critical phenomena, may be found in Ref. [19j. However, detailed 
studies of the strongly-coupled boson-fermion mixtures in the large A^ method are still missing. 

In this paper we present an extensive study of a large A' extension for a model of strongly-coupled 
boson-fermion mixtures originally proposed in Ref. [?]. We establish the l/A'-expansion in a theory 
of CFs which is equivalent to the original boson-fermion mixed system. We also derive an effective 
action of CFs up to the next-to-leading-order terms in the large A" expansion, and show that there 
can be the BCS superfluidity of CFs at low temperatures. 
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Our paper is organized as follows. In Sect. 2, we construct a large N extension of strongly-coupled 
boson-fermion mixtures at finite temperature and density based on the imaginary-time formalism. In 
Sect. 3, we rewrite the boson-fermion partition function in terms of CFs with the use of an auxiliary- 
field method. We derive an action functional of CFs and find a systematic expansion, l/A/'-expansion, 
which is equivalent to a loop expansion with respect to the CF fields. In Sect. 4, the l/A/'-expansion 
is employed to calculate the leading-order (LO) and the next-to-leading-order (NLO) terms in our 
CF action. We also derive a low-energy effective theory of CFs, and find that it reduces to a two- 
component free Fermi gas in the LO analysis, and to a weakly-interacting two-component Fermi gas 
up to the NLO study, which yields the superfiuidity of CFs at sufficiently low temperatures. Finally, 
in Sect. 5 we discuss the application of boson-fermiuon mixtures to dense QCD. In Appendix A, we 
give explicit forms of Fourier transformations especially for the proper vertex functions of CF fields. 
Appendix B provides details on the derivative expansion of the inverse propagator of CFs. 

2 Formulation of large N boson-fermion mixtures 

In our model, we treat bosons and two-component fermions using a nonrelativistic gas model of the 
boson-fermion mixture where bosons and fermions interact through a four-point contact interaction. 
We start from a Hamiltonian density of our boson-fermion mixture in three spatial dimensions, 

+ 1^ E E 4{^)^U^)M^)^^j{^) , (1) 

i,j = l cr=t,i 

where 0^ is the bosonic and tpcri is the fermionic field. We label the two internal states of the fermions 
by pseudospin indices a =t, | and extra large N indices of bosons and fermions by z, j = 1, 2, . . . , N. 
We assume that two different pseudospin states have the same mass and chemical potential (number 
density) and that the boson-fermion interaction is independent of pseudospin states. Setting = 1 
yields the same Hamiltonian density as in our previous work 7 , though here we neglect interactions 
between same species by assuming that the boson-fermion interaction is much stronger than the others. 
To make our analysis simple, we focus on an equally populated mixture of bosons and fermions, which 
means for each i we have n bosons and n fermions with an equal population in their number densities: 
Ubi = n^i + n^i = n and n^i = n^i = n/2. Also, we introduce total boson (fermion) number density 
as ntot = Nn. 

The bare boson-fermion coupling g^^ is related to the s-wave scattering length in the vacuum a^j 



3 



a, I 

— — >- 




(a) (b) 

Figure 1: Feynman diagrams for the boson- fermion mixture described by the action Eq.Q. The 
indices for large N degrees of freedom are exphcitly shown by i or j. (a) The single line denotes one 
of propagators for two-component fermions labeled by a and z, (b) the double line corresponds to the 
propagator for bosons, and (c) the interaction vertex between bosons and fermions is represented by 
the empty space associated with a suppression factor g^^jN in terms of the 1 /TV-expansion. 



by the following relation f20], 

1 f dk 1 

27ra^/ ~ g,f i|k|<A (27r)3 ^^(k) + £:/(k) ' 



(2) 



where Shi}^ = k /2m^ and S/(k) = k /2mj are the kinetic energies of the single boson and fermion, 
respectively, = m^m^/{m^ + mj) is the boson- fermion reduced mass, and A is a high- momentum 
cutoff of our model which sets a minimum atomic scale ro = (2A/7r)~^. For a simple notation, we 
will omit the constrain to the momentum-integral; |k| < A. 
The partition function at finite temperature becomes 

expressed by an imaginary-time functional integral over bosonic fields and fermionic Grassmann 

fields t/Jai^i^ai with the corresponding action functional of our boson- fermion mixture: 



i=l ^ 



i=l <T=t,4- ^ 



2m 



/ 



9. " 



Here for simple descriptions we have used notations: x = (x, r), Jdx = dr J dx, and h = ks = 1. 
Figure [l] shows the corresponding Feynman diagrams for the boson-fermion mixture, especially focused 
on the large N degrees of freedom. 

3 Prom boson-fermion mixtures to composite fermions 



In strongly-coupled boson-fermion mixtures characterized by the positive and small scattering length 
hf 



1 /3 

(0 < ^tot^bf ^ expect that bosons and fermions form bound dimers or composite fermions 
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(CFs) 4,5 and that low-energy phenomena can be described by an effective theory of these CFs. 
For this purpose, we introduce fermionic auxihary fields F^{x) and F^{x) by inserting the following 
identity into the partition function: 

1 = c /( n ^K^K) expj-^ [dxY, K{^)FUx)] , (5) 

where c is a normalization constant ^21^. We also define shifted fields F(j{x) and Fcr{x) as 



N 



^_bf_ 

N 



N 



^(l)i{x)llJai{x) + Fa{x) , 



N 



(6) 
(7) 



Note that the shifted fields F(j{x) and F(j{x) can be considered as fiuctuations of F^ and F^ around 
O(A^) singlet fields Xl^i ^i(^)'0cri(^)/^ and Ylf^i^*i{^)'4^(^i{^) respectively, both of which are 
arithmetic averages of many fields in terms of the large N . As we will see below we can in principle 
control these fiuctuations by changing N itself. The partition function Eq.Q then becomes 

Z = c lv[F'^,F'^mi^,]V[4^,i,^,,] 



exp< / dxdy 



N 



N 



^(j)1{x)D ^{x,y)(})i{y) ^ V^ai(^)^ ^{x,y)ilj^i{y) 

r- N 

L i,j = l cr=1-4 ^>'i <T=t4 ^ ^ 

L i=l i=l <T=t,^ 

r \ N ^ ^ 

+ dx — F„{x)F^{x) + Y,FAx)Mx)i^ai{x) + Y<PUxHai{x)F^{x) 

J — ^ I V •_i •_-] 



, (8) 



where D~^{x^y) and S~^{x^y) denote inverse Green's functions of bosons and fermions, respectively: 



D-\x,y) 
S-Hx,y) 



(9) 
(10) 



We also used a simple notation for functional integral measures: 

AT 



v[F,,F,mi^,m,,,^,,] = ( n ^f,^fM]Jv^;v^Mi[ H v^,,v^,X{ii) 
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Note that as in a usual Hubbard- Stratonovich transformation the introduction of auxiUary fields 
reduces the original bosonic and fermionic fields into bilinear forms which are diagonalized in terms 
of the large N indices, and they can be integrated out immediately. We first perform the fermionic 
functional integral in Eq.([8|, and the relevant part of the integration yields for each pair of indices 
<j ^ and i = 1, . . . , TV (i.e., summation over a and i is not assumed here), 

= j Vil)^iV^l)^ie^-pl^j dxdy[4)cji{x)^ j dzFcj{z)(l)i{z)S{z,x)\S~^{x,y)[^l)ai{y)^ j dwS{y,w)(l)*{w)F^{w)] 

— J dxdydzdwFcj{z)(j)i{z)S{z^ x)S~^{x, y)S{y^ w)(j)^{w)F(j{w)^ 
= j ^i^aPi^ai^^P^j dxdy i[j'^^{x)S~^{x,y)2lj'^^{y)'^ exp^^- j dzdwF^{z)(t)i{z)S{z,w)(l)1{w)F^{w) 
= exp|lndet[-5-i(x,^)]- J dxdy ^Ux)PAy)S{y.x)F,{x)My)'^ 
= exp|trln[-S'"^(x,?/)]-y" dxdy (l)*{x)A^{x,y)(l)i{y)^ , 

where "tr" and "det" are taken only over coordinate indices, and Acr{x^y) is defined by Acr{x^y) = 
F(j{y)S{y,x)F(j{x)^ whose ordering is important due to the anti-commuting nature of Grassmann 
fields F(j and F^-. We have also used a matrix formula: In det M = trlnM. Thus, Eq.|8| reduces to 



Z = c[zl{fif)] yp[F,,F,]P[(/>*,(/>,] 

exp<^ dxdy ^(jy'lix) D~^{x,y) - ^ A^{x,y) 
i=i L ^^^^^ 



(13) 



^{y) 



— jdxY, F,{x)F,{x)\ 



with a partition function for a two-component free Fermi gas Zl{ii^) = exp{2trln[—S'~^(x, ?/)]}. We 
proceed to perform bosonic functional integral, and the relevant part of the integration yields for each 
index i = 1, . . . , A^, 



exp< / dxdy (/>*(x) 



D \x,y)- A,{x,y) My)\ 
D-\x,y)^ ^ A,{x,y) \ 

det J dz[—D~^{x^ z)] S{z — y) — J dwD{z^w)A(j{w^y) 
cpj -tr In [-D~^{x, y)] -tr In 



(14) 



(12) 
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Applying Eq.(14) to Eq.(13), we obtain a partition function which is described only by CF fields: 



Z = c[4(mJ4(M/)]'^ 



(15) 



with a partition function for the ideal Bose gas ^0(^5) — exp{— tr ln[— I) ^(x, y)\ }. The correspond- 
ing action becomes 



S[F^,F^] = fdx Y F^{x)F^{x) ^ Ntiln S{x - y) - ^ j dwD{x,w)A^{w,y) 



.(16) 



Since N becomes an overall factor in the action and plays the same role as in a usual loop expansion, 
our 1/A^-expansion is equivalent to the loop expansion based on the CF action Eq.([l6]Q Let us 
normalize the CF fields as 

(17) 
(18) 
(19) 



T^{x) = VNF^ix), 
Fa{x) = VNF^{x), 
Aa{x,y) = Ta{y)S{y,x)F„{x) = NA^{x,y) 



which give an explicit form of the l/A'^-expansion, with the use of a formula in the logarithm: ln(l 

S[Fa,Fa] = fdx^ f'a{x)J'a{x)^NtT\n 5{x-y)-^ j dw ^ D{x,w)Aa{w,y) 

= — fdx ^cj{x)J^a{x) -Y(j^) T^^l Yl D{^^^)^CT 



{y)S{y,w)T^{w) 

(20) 



Here one can see that there is no internal degree of freedom associated with the large N extension 
and that 1/N only appears as a suppression factor of each higher-dimensional interaction between 
CFs. Then, we reach the following representation of the partition, 

Z = c'[Zl{^i,)Zl{^if)\'' jv[T„,Ta]ei,v{~S[Ta,Ta]} , (21) 

with a normalization constant c' . Figure |2] shows a formal expression of the CF action 5[^cr,^cr] of 



Eq.(20) in terms of Feynman graphs. We will give precise definitions of S ("self-energy"), r4 (4-point 



■"^In general when we have an overall factor 1/r in our action, propagators should be proportional to r, while any 
kind of vertices to 1/r. Then, any graph composed of P propagators and V vertices is proportional to r^~^ . On the 
other hand, such a graph has L = P — {V — 1) loops, which yields a relation: r^~^ = r^~^. Therefore, the series 
expansion in terms of r is equivalent to the loop expansion in diagrammatic expressions [22] , Setting 1/r = N yields 
our l/A/'-expansion. 
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Figure 2: Graphical representation of the CF action 5[^cr,^cr] in Eq.(20). The thick hnes represent 
the external lines of CF fields Ta and Ta- Definitions of S, r4 and V^n will be given later in Eq.(25), 
(48) and (67), respectively. 



vertex function) and (2n-point vertex function) later in Eq.(25), (48) and (67), respectively. It 
will be shown that the first two graphs in Fig |2] yields an inverse propagator of CFs which behaves as a 
free Fermi particle within our approximation, and the rest of graphs can be considered as interaction 
vertices of CFs. 

4 1/A^-expansion of strongly-coupled boson- fermion mixtures 



In the following, we will perform the 1/A^-expansion based on Eq.(20) up to the next-to- leading- 
order terms and derive a low-energy effective theory of CFs. We will show that under an assump- 
tion discussed below an effective interaction between CFs are weakly attractive, and that the BCS- 
superfiuidity of CFs (CF-BCS) is realized at sufficiently low temperatures. 

4.1 The leading-order terms 



The leading-order (LO) terms in the 1/A^-expansion, i.e., 0(1) terms in Eq.(20), become 



^'Lo\^a-,^a\ = — — f dx ^ Jv(x)J>(x)— j dxdw ^ J^^{x)D{x^w)S{x^w)J^^{w) ^ (22) 



which we can rewrite with the use of Fourier transforms (see Appendix [A|) as, 

9 u f J ^ ^ I J ^ ^ \ - J 



J'Ap) 



= - jdpY. ^Ap)G-\p)T,{p) . (23) 

Here we have used notations: p = (p, iuj)^ Jdp = T uj Jdp/ (27r)^ with the Matsubara frequency uj 
and spatial momentum vector p. Also, we have introduced an inverse propagator of CF fields J> as 

G-\p) = — + S(p) , (24) 

with a CF "self-energy", or single "bubble" of bosons and fermions, S(p) given by 

S(p) = JdqD{q)S{p-q). (25) 



(J^ i 
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Figure 3: Graphical representation of the propagator for "bare" CF fields F(j (not for Jv). In the 
right hand side, the external triple line denotes a pair of single boson and fermion, and only appears 
as a shorthand notation for a half of the vertex in Fig.[l|c). As for the propagator of J^a{'= VNE^j) 
fields, we need to multiply both sides by the normalization factor N . 



Equations ( 24 ) and ( 25 ) show that the propagator is represented by an infinite geometric series of the 



original boson- fermion bubbles, as shown in Fig.jsj In the right hand side of Fig.|3| we can see that the 
n-th graph has large N power-counting factors (i) N'^~'^ from n — 1 internal loops, (ii) (l/N)'^ from n 
vertices, and (iii) (a/TV)^ from the normalization (J> = ^/NF^j) for any n G N, which give an 0{1/N^) 
term in total, i.e., the leading-order contribution in the l/A/'-expansion as shown in Eqs.([23|)-([25|). 
Let us expand the inverse propagator in order to derive a low-energy effective theory of CFs. 



The summation over the bosonic Matsubara frequency in the self-energy Eq.(|25| is performed as 



j dqD{q)Sip - q) 

f dq 1 



1 



-I 



(27r)3 iuj^ - ^^(q) i{uj - uj^) - ^^(p - q) 

e^^ 1 1 



dq 1 
lim 



dz- 



(27r)3 ^^lo 2711 Jc ef^^ - 1 z - ^(q) icj - z - ^^(p - q) 
dq 1 - n/(p - q) + n5(q) 



(26) 



(2^)3 e/p-q)+4(q)-zc^' 
where n5(k) and n/(k) denote the Bose-Einstein and Fermi-Dirac distribution functions, respectively, 

1 



715 (k) 

nf{k) = 



1 



(27) 
(28) 



with and kinetic energies of single boson and fermion relative to the chemical potentials /i^ and 
/ij, respectively: ^^(k) = k^/2m^ — /i^ , ^j(k) = k^/2mj — Mj- Here we have taken a standard contour 



Imz 




Re 2; 



Figure 4: Contour for evaluation of the Matsubara frequency sum in Eq.(26). The crosses indicate 
the position of the poles in the integrand. 

C on the complex z-plane (see Fig.[4| in order to convert the summation over the bosonic Matsubara 
frequency cj^ into a complex-integration along C 22 . Then the inverse propagator Eq.(24) reads 

1 1 - n/(p - q) +n5(q) 1 



G-\p) 



27ra 



W U(q) 6(p-q) + 4(q)-^^ 



(29) 



where the kinetic energy in the relative coordinate £:(q) is defined by £(q) = q^/(2m^), and by using 



Eq.(|2| we replaced the coupling constant g^^ with the scattering length a^j. Note that up to this 
stage there is no need to put any assumption related to the strength of our coupling constant g^^ , or 



1/3 



1 /3 

Now let us study a strongly-coupled mixture characterized by < n^'^^ a^^ <C 1. In this case, it 
is natural to consider the situation that both /i^ and /ij are almost equal to — a;^^/2, where uo^^ is 
a binding energy of an isolated boson- fermion pair in the vacuum: uj^^ = l/(2m^a^j). This is in 
accordance with the fact that in the strongly-coupled mixture the system becomes a dilute gas of CFs 

1 /3 

due to n^Q^ \%f \ <^ 1- We will later see that the number equations to relate the chemical potentials 
and the particle density indeed have a solution /i^ + /i^ ~^bf ' "^^^^ implies putting one more pair 
reduces the total energy by an energy almost equal to a; . 



Then, the low energy and low momentum expansions of Eq.(29) at zero temperature gives (see 
Appendix p| 



47r 



(30) 



%f V 2(m^ + m^) 

where ji denotes a total boson- fermion chemical potential: /i = (i^ ^ ii^{< 0), and the chemical 
potentials are yet to be determined. We proceed to expand Eq.(30) in terms of {p^/[2(m^ + ^/)] ~ 
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to obtain the derivative expansion of the inverse propagator, 

G-\p,E) a-cp^+d^, (31) 
with the zero temperature coefficients 



27ra^j 27r 



1 f^R 

47r rrij + \j 2|/i 



(32) 
(33) 



d = 7:A- (34) 



(35) 



Thus in low-energy scales, Eq.(23) can be approximated by its effective action, 

^LO i^a^J^a] ^ - Yl ^^^(P) (^^^ - 2 [d/{2c)] ^ ^^<^ijP) ' 

Performing a proper normalization of CF fields with = VdJ^a- , = \fdTa yields a low-energy 
effective action, 

sfo = - E (^^ - ^ + ^-^^'^ ' (=^6) 

with the kinetic mass = d/{2c) and the chemical potential /ip = a/dioi the normalized CF fields 



mp = mj+m^, (37) 
= 2c^^^(vT-X). (38) 



Here we have defined a dimensionless parameter X = \/i\/uj^_^. Using the effective action Eq.(36), we 



can construct an effective theory described by the following partition function, 

Zfo = c[Z'MZf{^f)f ff]J V^^V^X (39) 

which is valid for phenomena dominated by low-energy and low- momentum scales such that uj/\/j.\ <C 
1, £:(p)/|/i| <C 1, and the chemical potential is determined by number equations, 

ntot = -^InZ = --^InZ. (40) 

We will estimate chemical potentials, /ij, /i^ and /ip, at zero temperature, where the number 
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0.25 0.50 0.75 1 

Figure 5: Numerical plot of /(X) = (VX - X)(2 - \|^fXfl^ as a function of X. 



equations reduce to 



^tot 



d 
d 



InZ 



InZ; 



■eff. 
LO 



d 



37r2 



■eff. 
LO 



1 

^X 



-X 



dk 

(2^ 



■n5(k) 



(41) 



Here we used the fact that n5(k) vanishes at T = with ja^ < and also that the effective action 



Eq.(36) is the same action as for a two-component free Fermi gas with a mass and a chemical 



potential /ip. We can rewrite Eq.(41) as a dimensionless equation 



1 



=■■ fix), (42) 

where we defined a function / as /(X) = (a/X — X)(2 — 1/a/X)^^^. Figure [s] shows a numerical plot 
of /(X) as a function of X, and we can see that /(X) becomes zero at X = 1/4 and X = 1. Since 



the left hand side of Eq.(42) becomes quite small in the strongly-coupled mixture, Eq.(42) will give 



two solutions around X ~ 1/4 and X ~ 1. From now on, we will focus on a solution X ~ 1 which is 
consistent with the case of the dilute gas as we have mentioned before. 

Based on the above argument, we introduce a small positive parameter S{= 1 — X) which should 



be determined through the number equations. Putting this into Eq.(|41|) yields 

(STr^ntot)'/^ 



2m, 



= X 2 



2/3 



(43) 
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which finahy gives 



(44) 
(45) 



with a Fermi energy of CFs ep = (37r^ntot)^^^/(2mp). Note that the above analysis becomes rehable 
^ cf/uj^^) whi( 



only with a small ep/uj^^) which demands the following condition: 



2m, 



(37r^ 



2x2/3 \/ 1/3 



(46) 



Now we can see that Eq.(46) is automaticahy satisfied in our strongly-coupled mixture such that 



< n\ila,j « 1. 



Combining Eqs.(36) and (45), we find that in the leading order of the l/A/'-expansion a low- 



energy effective theory of the strongly-coupled boson-fermion mixture, whose energy scales satisfy 
"^/^bf ^ 1' ^l^hf ^ ^ ^{p)/^t,f ^ 1? becomes just a two-component free Fermi gas of CFs 
with a mass = + and the same number density as the total number density of the original 
fermions ntot- We remark that in our strongly-coupled mixture the small expansion parameters are 
1/N and n^of a^^. 

4.2 The next-to-leading-order term 



We proceed to study the next-to- leading-order (NLO) term in Eq.(20), and we will show that the 



NLO term gives an effective attraction between CFs in low-energy scales. The NLO term in Eq.(20) 
becomes 

= f dvdwdz ^ D{x,w)Aa{w,v)D{v^ z)Ap{z^y) 

= ^ I dxdvdwdz D{x^w)Fa{v)S{v^w)Fa{w)D{v^z)iFp{x)S{x^z)Fp{z) 

X j dqD{q)S{p2 - q)D{pi - p2 -^q)S{p3 - q) ^a{Pi)^p{P3)J^p{P4)J^a{P2) 

^ \ "Y [{ n^^i)^(^4 +P2 -Pl -Ps) ^A{{Pi}Ul)^cT{Pl)^p{P3)J^p{P4)J^c 



iP2), (47) 



where r4 represents the proper 4-point vertex of CFs, defined by 

r4 (telti) = -^Jdq D{q)S{p2 - q)D{p^ - P2 + q)S{ps - q) , 



(48) 
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P4;p P3;p P4;p ps-q]p,i Ps^p 

Figure 6: Graphical representation of the proper 4-point vertex of "bare" CF fields Fcr (again, not for 
J-a)' In order to obtain the proper 4-point vertex of Jv, we need to multiply both sides by (a/TV)^, 
which comes from the normalization of four external legs (J> = ^fNFa). 



with a set of momenta {pi}^^i = {pi, Ps, Pa} constrained by the energy-momentum conservation, 
P4 -\- P2 — Pi — Ps =0. The graphical representation of Eq.(48) is shown in Fig.[6j which indicates that 



the CFs interact through the exchange of their constituent particles. We can see that the right hand 
side of Fig.[6]gives large N power-counting factors (i) N from an internal loop, (ii) from four 

vertices, and (iii) (a/TV)^ from the normalization (J> = y/NFa-)^ which yield an 0{1/N) term as is 
r4 in Eq.(|48|. 

For a low-energy effective theory, we expand the proper 4-point vertex in terms of E /uj^^ with a 
small energy scale E relative to its binding energy uo^^ . Then the dominant contribution becomes 

^i{{Pi}) = r4({(0,7rT),(0,7rT),(0,-7rT),(0,-7rT)}) + 0(i?/w^^). (49) 

We note that at finite temperatures it is impossible to put the frequencies in r4 equal to zero, since 
the CF fields in Eq.(|47| are Grassmann fields (purely fermionic) and do not have Matsubara zero 
mode, i.e., uj = (2n + l)7rT, n G Z. Note also that our procedure is essentially the same as in the 
derivation of the effective theory for Cooper pairs, which can be considered as composite bosons in 



two-component Fermi gases 16-18,23 26 . Let us denote the low-energy effective vertex in Eq.(49) 
by r4(0), which becomes 

r4(0) = r4({(0,^T),(0,^T),(0,-^T),(0,-^T)}) 

= ~N Jj^^T. - iw, - C/-q) -zttT - iw, - C/-q) (iw, - 4(q)) ' 
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Imz 




Rez 



Figure 7: Contour for evaluation of the Matsubara frequency sum in Eq.(51). The crosses indicate 
the position of the poles in the integrand. 



We can perform the summation over the bosonic Matsubara frequency uj^ in Eq.(50) as follows, 

2 



1 



1 



1 



d 



2m Jc e^^ - 1 z-iTiT^ ^^(q) z + zttT + ^^(q) ds \z - ^^(q) - s 



1 



= [1 -n/(q) +n5(q)]- 



2[^/q)+4(q)] 



1 



=0 
(q) 



. (51) 



{[^/(q) +4(q)]' + (^T)2}' [^/q) + 4(q)]' + (^^)' ^^^^^^ 
Here as in the calculation of Eq.([26|, we have taken a standard contour C on the complex z-plane (see 
Fig.[7|) in order to convert the Matsubara summation into a complex- integration along C. Substituting 



Eq.([51| into ([50|) yields 

r4(o) 



1 f 2[^/q)+4(q)] [l- nj(q) + n5(q) 

2 , / rx.XOl 2 



TV 



h^^f {[e,(q) + 4(qr+W}^ 

^ f_M 1 5nb(q) 

( 



(52) 



A^y(27r)3[^^(q)+4(q)]2+(^r)2 5e,(q) " 

We will estimate r4(0) analytically at T = with the same assumption as before, that is, with an 
assumption that both /i^ and /ij are negative and their magnitudes are almost equal to ^f,^/2. This 



assumption still holds since the NLO term Eq.(47) is suppressed by a factor 1/A^ compared to the LO 



terms and does not change the values of chemical potentials so much from Eq.(44) and (45). Under 
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this assumption, r4(0) at T = reduces to 

r = -1 f ^ 



1 (2m^)3/2 /-A/v^ 



Jo 



N 7r2 Jo (x2 + |/x|)3 

l(2mj3/2ir 1 / A 



dx 



(A/V2^)^ - |M|A/(2mJ 



Ml/ H[|M|+A2/(2mJ]^ 



1 1 / m 



A/'27r V2|/i 



3/2 



(53) 



where we have neglected 0{y^2m^ l/^l/^) corrections in the finial step. Thus in low-energy scales, 



Eq.(47) can be approximated by 

5nLO [^a^ J^a] 



(P2) . (54) 



According to the analysis on the LO terms, we perform the same normalization as in Eq.(36), = 
\fdTa^ ^cr = \fdTa^ which yields a low-energy effective action in the NLO, 

= \ ^ /'(n#i)<^(P4+P2-pi-p3)5pp*a(pi)tp(p3)*p(p4)*a(P2). (55) 



Here we defined an effective four- Fermi coupling constant by g^^ = r4{0)/(P, which becomes 

_ 1 27rafc^ 



(56) 



We apply the LO result Eq.(44), i.e., X ^ 1 — e^/cj^^, to the above Eq.(56), and finally obtain 



1 27Ta^j 
~ TV * 

Introducing the density of states per unit volume at the Fermi surface Nf{0): 

we have the dimensionless parameter for the strength of the effective four-Fermi interaction. 



(57) 



(58) 



NF{0)g, 



= -2 



1/3 ^1/3^ 

7T J N 



(59) 
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Here we find that the NLO term yields an effective four-Fermi interaction, which is attractive and 
weak in a twofold meaning. First, we are considering strongly-coupled boson- fermion mixtures so 

1 /3 

that the dimensionless parameter n^^^ a^^ is positive and much smaller than 1, which makes NF{{))g^^ 
negative and much smaller than 1, as discussed in Ref. 7 ioi N = 1. Secondly, we also have large 
N degrees of freedom in the original boson- fermion mixture, which yield the large suppression factor 



l/N in Eq.(59). 



4.3 BCS superfluidity of composite fermions 

From the results of the previous sections, we have a low-energy effective action up to the NLO term 
in the l/A/'-expansion: 



dp crip) [ ioj 



2m^ 



I ^ [ ( n^O'^^^^ +P2 -Pi -ps) ^pp^a(Pl)^p(P3)^p(P4)^a(P2) 



(60) 



with the physical parameters of CF fields. 



Mf 
9ff 



2m^ 



1 27ra 



bf 



N 



(61) 
(62) 

(63) 



Equation (60) is nothing but an action of two-component Fermi gases with weakly attractive four- 



Fermi interactions, which yields the BCS-paired state of fermions at low temperature. Thus we can 



expect that our system described by Eq.(60) also favor the BCS superfluidity of CFs (CF-BCS) below 
a transition temperature 



27,28 



Tc (CF-BCS) 



- - e^exp 
TT \ e J \ zhpapp 



(64) 



with the Fermi momentum = yjlm^ep = (37r^ntot)^^^ and the s-wave scattering length given by 



777^ 

47r' 



1 

"7V2m/^^ 



(65) 



Setting = 1 yields the same result as derived in our previous work 7 . Note that the coefficient of 
a^j in Eq.(65) becomes —2 for A/" = 1 with — m^, which is the same in magnitude but opposite in 
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Figure 8: Graphical representation of the proper 6-point vertex of "bare" CF fields F^. In order to 
obtain the proper 6-point vertex of Jv, we need to multiply both sides by (a/TV)^, which comes from 
the normalization of six external legs (Jv = VNFct). 

sign from the scattering length between bosonic dimers composed of spin-singlet fermion pairs within 
the same approximation. This is because our CFs are different in the statistics of their constituent 
particles from the composite bosons, so-called Cooper pairs, in two-component Fermi gases (5[ [23|[24| . 

4.4 Higher order terms in the 1/A^-expansion 



We now consider the higher order terms in the 1/A^-expansion of the CF action Eq.(20), especially 
focusing on the corrections to the effective four-Fermi. 

Figure [8] shows a graphical representation of the proper 6-point vertex, based on the sextet term of 
CF fields in Eq.(20). We can see that the right hand side of Fig. [s] gives large TV power-counting factors 



(i) N from an internal loop, (ii) from six vertices, and (iii) (vA^)^ from the normalization 



= \/NF^), which yield an 0(l/N'^) term as it should be in Eq.(20). 



In general, we denote an 0{1/N'^ ^) term in Eq.(20) by S^"^ (n > 2), which is composed of 2n 
CF fields and the proper 2n-point vertex, 

n-l 



1 f f ^ \ f ^ ^ 

= - / ( n *2,-i^iP2, ) ^P2, - 



X r2n(te}£l)| n/'.,(p2,-l)^.,(p2,)| . (66) 



Here r2n(te}i=i) represents a proper 2n-point vertex of CF fields, defined by 

^2n{{p^}tl) = dql[Di ^P2.-l -^to+g S( ^^2.-^^-1-^ ,(67) 

j = l ^ /c=l ^ ^ k=l k=l ^ 

with a set of momenta {j>i\\^\ — {PI7P27 • • • ^V2n\ constrained by the energy- momentum conservation, 
Xlj=iP2j — Xlj=i^'2j-i = 0- Appendix [a| gives a detail derivation of the above expressions. Figure 
[9] shows a graphical representation of the proper 2n-point vertex. We can see that the right hand 
side of Fig. [9] gives large N power-counting factors (i) N from an internal loop, (ii) (XjN^^^ from 2n 
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Figure 9: Graphical representation of the proper 2n-points vertex of "bare" CFs F^r. In order to 
obtain the proper 2n-points vertex of Jv, we need to multiply both sides by (V^)^"^, which comes 
from the normalization of 2n external legs (Jv = VNFa)- 



vertices, and (iii) [yNY^ from the normalization (Jv = vNFfj)^ which yield an 0(1 /N^ ^) term as 



is in Eq.(p7|). Comparing Eqs.(66), (67) with Eqs.(47), (48), we can easily check the consistency 
for the case of n = 2, i.e., S^^^ \_?(t-,^(t\ = 5nlo [-^a, -^a] • 

Now let us consider diagrams higher order in the l/A/'-expansion which contribute to the effective 
four-Fermi interaction. The leading-order contribution to the effective four-Fermi interaction is given 
by the proper 4-point vertex Since our l/A/'-expansion is equivalent to the loop expansion in terms 
of CF fields, the corrections to the proper 4-point vertex T/^ should start from one-loop graphs, as 
shown in Fig.p^Q] and pTj Figure [lO] represents the contraction of the proper 6-point vertex Fg, and 
its right hand side gives large N power-counting factors (i) N from an internal loop, (ii) (1/N)^ from 
six vertices, and (iii) {s/N)^ from the normalization (J> = ^/NF^j)^ which give an 0{1/N'^) term in 
total, i.e., the next-to- leading-order contribution to the effective four-Fermi interaction. Figure 11 
shows the shortest ladder diagram of the proper 4-point vertices and its right hand side gives 
large N power-counting factors (i) N'^ from two internal loops, (ii) from eight vertices, and 

(iii) [s/N)^ from the normalization (Jv = \/NFcj)^ which again yield an 0{1/N'^) term. 

It is also possible to consider general 0(1 /N^) corrections based on the loop expansion in terms 
of CF fields which contribute to the effective four- Fermi interaction. Instead, we just remark that all 
the graphs are in the same order in terms of n^/^a^j, thus for A^" = 1 we need to sum up them to 



obtain the effective vertex function, as performed numerically in Refs 29 - 31 . 





Figure 10: Graphical representation of the contraction of the proper 6-point vertex of "bare" CF fields 
Fcr, which contributes to the effective four-Fermi interaction. 
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Figure 11: Graphical representation of the shortest ladder diagram of the proper 4-point vertices of 
"bare" CF fields F^r, which contributes to the effective four- Fermi interaction. 

5 Summary and discussion 

We have investigated the large N expansion for strongly-coupled boson-fermion mixtures, proposed in 
Ref. [7]. We first derived a theory equivalent to the original boson-fermion mixture, which is described 
by composite fermions (CFs). The 1/A^-expansion naturally appears in the quantum theory of CFs. 
We showed that the leading-order terms in the 1/A^-expansion yield a low energy effective action of 
CFs which is equivalent to that of a two-component free Fermi gas. The next-to- leading-order term 
was also estimated, and it turned out that the effective action up to the NLO reduces to an action of a 
weakly-interacting two-component Fermi gas. Thus we concluded that there is the BCS superfiuidity 



of CFs below Tc(CF-BCS) given by Eq.(64) in our large N model. Also we discussed how to estimate 
the higher order terms in the 1 /TV-expansion, where the diagrammatic representation provides simple 
explanations for power-counting of the 1/A^ factors. 

Finally, we would like to mention important similarities between our boson-fermion mixtures and 
hadron physics. Table [l] summarizes the correspondence in components between ultracold atoms and 
dense QCD, both of which can be considered as boson-fermion mixtures with = 1 in our model [6]-[8]. 
Note that it is known that the effective interaction between nucleons is not so strong as the original 
gluonic interaction between quarks. For example, the energy gap in nuclear matter, i.e., superfluid 
matter of nucleons, is at most a few MeV 32 , while the energy gap in color-superconductivity is from 
10 to 10^ times larger 33 . This is consistent with our results which show that the weakly-coupled 
CF system can be derived from the strongly-coupled boson-fermion mixture. However, our model 
seems too simple to relate its results to various phenomena in QCD. Furthermore, chiral symmetry 
breaking plays an important role in hadron physics 6 , which does not appear in nonrelativistic 
systems. Keeping these observations in mind, we suggest that both theoretical and experimental 
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Table 1: Correspondence between the boson- fermion mixture in ultracold atoms and the diquark- 
quark mixture in dense QCD. 



Our Notations 


Ultra- Cold Atoms 


Dense QCD 




bosonic atom (e.g.,^^Rb) 


diquark 




fermionic atoms (e.g.,^^K) 


unpaired quarks 




composite fermions (boson-fermion dimers) 


nucleons (baryons) 




boson-fermion attraction 


gluonic attraction 




composite-fermion superfluidity 


nucleon superfluidity 



studies in boson-fermion mixtures provide a new tool to investigate properties of dense QCD, which 
is not readily observable in laboratory experiments. 

Acknowledgements 

The author thanks G. Baym, D. Blaschke, T. Hatsuda, S. Uchino for fruitful discussions. This research 
was supported by JSPS Research Fellowship for young scientists. 

A Fourier transformations 

The definitions of Fourier transforms used in our main text are given by 

Ta{x) = jdve'^P^T^ip) , (68) 



j;(a;) = J dpe-'P^T^{p), (69) 
D{x,y) = Jdqe''''-''-y'>D{q) , (70) 
S{x,y) = [dqe"^^--y^S{q), (71) 



with an inner product: ipx = iur — ip • x where uj and p denote the Matsubara frequency and 
spatial momentum vector, respectively. For simplicity, we have used notations: p = (p,ia;), Jdp = 
T^^uj Jdp/{27r)^. We adopt a convention which distinguishes functions from their Fourier trans- 
forms only by their arguments. Correspondingly, the Fourier transform of the quadratic term in 



Eq.(|22| becomes 

dxdw J'a{x)D{x^ w)S{x^ w)T(j{u)) 

J dpidp2dqidq2 S{qi q2 - Pi)S{p2 - qi - q2)f'a{Pi)D{qi)S{q2)Ta{P2) 



jdp^^ip) JdqD{q)S{p-q) 



Mp) ■ (72) 
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Also, the quartic term in Eq.(47) transforms as 



J dxdvdwdzD{x^w)J'cr{'^)S{v,w)J'cr{uj)D{v,z)Tp{x)S{x,z)J^p{z) 

J (^dpidqi^S{qi +^4-^3)^(^2 - qi - q2)S{q2 ^ Qs - Pi)S{pa - qs - q^) 

= j []\dV^5{VA^P2 -Pl -P3)^{{Pi}Ui)^a{Pl)^p{P3)J^p{P4.)J^a{P2) , 

where r({p^}|^^) represents a proper 4-point vertex of composite fermions, defined by 
r({Pi}tl) = JdqD{q)S{p2-q)D{pi-p2 + q)S{ps-q), 

with a set of momenta {pi}j^i = {^1,^2,^3,^4}- 

We can write down a general 2n-point vertex function explicitly in its Fourier transform, 

1 



(73) 



(74) 



TV" 



-tr 



Jdw^ D{x,w)Aa{w,y) 



(75) 



^ ^ ^ 2n+l N ^ n 

— n ^ / ( n - ^2n+l)s YlD{x2j-l,X2j)f'a,{x2j + l)S{x2j^l,X2j)Taj{x2j) 



Y[ dq2j-idq2jdp2j-idp2j 



S{qi^q2n-p2n-i)l Yi S{p2j - q2j - q2j -i) \ I Y{S{q2j^i^q2j -P2j-i)\ 
I j=i J ^ j=i J 

x(f[D{q2j-i)S{q2j)\\ n^-.(^2,-i)^.,(p2,)| 
I j=i ^ ^ j=i ^ 

S / X]^2/c-l - ^P2k^qi ] si ^P2k - ^P2k-1 -qij >< n'^^^(^2j-l)^a,(P2j) [ 

( I[dp2j-ldp2j] S( ^P2, -^P2,-l)r2n(te}-"i)<^ H ^^-^ (^2, - 1 ) ^a, (P2, ) [ , 



where ^2n{{Pi}i=i) represents the proper 2n-point vertex of composite fermions, defined by 

r2n(te}-^l) = -^^T^iri / ^^n^( - ^P2k^q \ si ^P2k - ^P2k-l - ^j , 

with a set of momenta {pi}1^i — {pi,P2, • • • iP2n\' 



(76) 
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B Derivative expansion of the inverse propagator for compos- 
ite fermions at T = 



We will give details on the derivation of Eq.(30) and Eqs.(32)-(34) in the derivative expansion of 



G ^{p) at zero temperature. We can formally perform a real-time analysis by replacing our Matsubara 
frequency iuo with a continuous energy variable E at zero temperature. Using a four-momentum in 



real-time formalism; (p,£^), we have a real-time form of Eq.(29) at zero temperature, 

1 





f dq_ 




J (2^)3 




f dq 




J {2nf 


^R 





1 



1 1 



^(q) ^(q) - p • qM/ + pV(2^/) - ^ - + /^/) 



(77) 



Here we used the fact that n5(q) and n/(p — q) vanishes at T = under our assumption: /i^ < and 



/ij < 0, and we denoted an integral in Eq.(77) by /i, which becomes 



J (2^; 



1 



q dq 
1 

4^ 



1 1 

^(q) ^(q) - p • qM/ +^7(2^/) - e 
1 1 



1 I ^(q) e[ci) — pq cos O/rrij^ -\- A 



q^dq 

pq 



In 



-pq/rUf + e{q) + A 



(78) 



with p = |p|, q = |q|, /j, = fif{< 0), and A = p^ / (2m — E -\- \fi\. As we will see below, Ii yields 



a finite value even in the limit of A ^ oo. We can rewrite terms in the logarithm in Eq.(78) as 

^2 



PQ 



+ £:(q) + A = 



1 



2m^ 
1 

2m ^ 



2 o^i? 

rrir 



pq 



p- 



2m 



f 



{{q^apf^2m^A'}, 



and 



PQ 

mr 



+ £:(q) + A = 



2m, 



-{{q-apf^2m^A'}, 



(79) 



(80) 



with a mass-ratio parameter a = m^/mn^ and A' =p^ /[2{m^ m r)] — E ^ Then, Eq.(78) reads 



h 



m 



m^A 



/ 



dq q In 



{q + ap)'^ + 2m^A' 



{q - apY + 2mj^A' 



By employing the following integral formula: 



dx X In 



(x - a)2 + 6 



(81) 



(82) 



In 



(x - a)2 + 6 



2ax — 2ay\b\ < tan 



- tan 



X — a\ 1 
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we can perform the integration in Eq.(l78|), 



In 



(A + ap)'^ + 2m^A' 



Att'^P \ 

l^2apA + ^^2apV|2m^A'||tan-^ 



(A - + 2m^A^ 
A-\- ap 





j 










27r2 


27r2 



In 



(A + ap) 



tan 



A — 



^)} 



2 ^ 2m^A' 



(A - + 2m^A' 



f V|2m^A^||tan-^ 



A + 



- tan 



A — ap 



(83) 



Noting that the cutoff A is the largest scale in our model and also that we are interested in the low- 



energy and low- momentum regime, let us expand the logarithm and arctangent functions in Eq.(83). 
As for the logarithm function, the power-series formula, ln(l -\- x) = l)"^"^^;"^/^, gives 



A^ - aV + 2m^A' 



In 



(A + ap)^ + 2m^A' 



(A - ap)2 + 2m^A^ 



A2 
2 
A 



a^p^ - 2mj^A' 
A2 

a^p^ — 2m^A' 
A2 



In 



1 + 2ap/A + (aV + 2m^A')/A2 



1 - 2ap/A + (aV + 2m^A0/A2 
ap{a'^p'^ + 2mj^A') 



Aap^O 



A2 



= 2apA + 0("^("^^^f^---^^)), 

while the formula for the arctangent function: tan~^(x) = 7r/2 — tan~^(l/x), yields 

1 / A + ap \ _i / A — ap 

' ' ' tan ' 



tan 



which finally give an explicit form of Ii within the derivative expansion, 

' ap{a^p^ ± 2m p. A') 



rrij^A rrij^A rn^A 
^ ^ ~ ~2^ ~ ~2^ 



27r 



n th 

\/|2m,A 



O 



A±ap 



(84) 



(85) 



O 



A 



O 



A ± ap 



(2mj3/2 



47r 



\A' 



O 



ap{a p ± 2m^A') 



A 



A ± ap 



Substituting Eq.(86) into Eq.(77) yields 



27ra 



47r 



2(m^ + m^) 

.2 . 



(87) 



which finally gives the derivative expansion of G ^ in terms of {p^/[2(m^ + ^/)] ~ ^}/ImI under the 
assumption sated in the main text. 



(2m^)3/2 



27ra 



bf 



m 



27ia 



^ 1 
1 ^ 



a 



bf 



pV[2(r 



1 



27r2(m^ + m^) Y 2|/i| 



27r V 2|/i| 



The above expression corresponds to Eq.(31) with the coefficients a, c, (i given by Eq.(32)-(34). 
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